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ABSTRACT 

The spike (S) protein of corona viruses, a type I membrane glycoprotein, is primarily responsible for entry into 
susceptible cells by binding with specific receptors on cells and mediating subsequent virus-cell fusion. The 
bovine coronavirus (BCoV) S protein is cleaved into two subunits, the N-terminal SI and the C-terminal S2. 
The proteolytic cleavage site of S protein is highly conserved among BCoV strains and is located between 
amino acids 763 and 768 (KRRSRR). This study describes a single mutation in the S protein cleavage site of 
three Brazilian strains of BCoV detected in diarrheic fecal samples from calves naturally infected. The sequenced 
PCR products revealed that amino acid sequence of the cleavage site of our strains was KRRSSR, indicating 
a mutation at amino acid position 767 (R—> S). This amino acid substitution occurred due to a single nucleotide 
substitution in the sequence of DNA corresponding to the proteolytic cleavage site, CGT to AGT. This is the 
first description of this nucleotide mutation (C to A), which resulted in the substitution of arginine to serine 
in the S cleavage site. In this study we speculated the probable effects of this mutation in the proteolytic 
cleavage site using the murine hepatitis coronavirus (MHV) as a comparative model. 
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Bovine coronavirus (BCoV), a member of the family 
Coronaviridae, order Nidovirales, belongs to group 2 of the 
coronaviruses which include murine hepatitis coronaviruses 
(MHV), porcine hemagglutinating encephalomyelitis virus 
(HEV), equine coronavirus (ECoV), rat coronavirus (RtCoV) and 
human respiratory coronavirus (HCoV-OC43) (9). BCoV is an 
enveloped virus with single-stranded, positive-sense RNA 
genome of approximately 32 kb length that encodes five major 
structural proteins: the nucleocapsid (N), the transmembrane 
(M), the hemaglutinin esterase (HE), the spike (S), and the small 
protein (E) (13). The S glycoprotein of BCoV is a large membrane 
glycoprotein of approximately 150 kDa that forms the peplomers 
(club-shaped structures) on the virion surface. The S protein is 
primarily responsible for the entry of coronavirus into 
susceptible cells by binding to specific receptors on cells and 
mediating virus-cell fusion and subsequent cell-cell fusion 
during infection (5). 


In several coronaviruses, such as infectious bronchitis virus 
(IB V), MHV and BCoV, as a late event in maturation, the protein 
is cleaved into two subunits: S1 (aminoterminal region) and S2 
(carboxyterminal region) (17). Proteolytic cleavage of the S 
protein of these coronaviruses occurs adjacent to a sequence 
of basic amino acids on the carboxyterminal region of S1. In the 
S protein of BCoV, a predicted basic amino acid sequence 
(KRRSRR) is involved in the cleavage by the host cell-derived 
proteolytic enzyme. This sequence, highly conserved among 
BCoV strains, encompasses amino acids 763 to 768; the cleavage 
occurs between amino acids 768 and 769 (19). The cleavage of 
the S protein has been reported as a process related to the viral 
infectivity and cell fusion from other group 2 coronaviruses. 
Studies related to MHV, the best-studied member of the 
Coronavirus family, demonstrated that cleavage of S is not 
essential for infectivity but is associated with enhanced cell 
fusion (syncytia) in infected cell monolayers (7,8,20). 
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This study describes a single mutation in the S protein 
cleavage site of wild type Brazilian strains of BCoV detected in 
calves naturally infected and speculates the possible biological 
effects of this mutation at the proteolytic cleavage site using 
the MHV as a comparative model. 

Three BCoV positive fecal samples (BR-UEL1, BR-UEL2 and 
BR-UEL3) were obtained from calves up to 30 days old with 
clinical signs of diarrhea in a Brazilian dairy herd from Minas 
Gerais State (21°4T49” S; 45° 18'45” W). These samples were 
previously identified as BCoV positive by RT-PCR assay for N 
gene detection (21) and negative for bovine group A rotavirus 
and Cryptosporidium sp by polyacrylamide gel electrophoresis 
technique and modified Ziehl-Nielsen method. Fecal samples 
were diluted 2-fold in 0.01 M phosphate-buffered saline (PBS) 
pH 7.2 (137 mM NaCl; 3 mM KC1; 8 mM Na 2 HP0 4 ; 15 mM 
KH 2 HP0 4 ), centrifuged at 3000 x g for 15 min at 4°C to remove the 
cell debris. The supernatants were used for RNA extraction. 
Aliquots of 400 pi from fecal suspensions were treated with SDS 
at a final concentration of 1 % (w/v), homogenized by vortexing 
and kept at 56°C for 30 min. A combination of phenol/chloroform/ 
isoamyl alcohol and silica/guanidinium isothiocyanate methods 
was performed according to Barreiros et al. (1) with slight 
modifications. Briefly, 400 pi of phenol/chloroform/isoamyl 
alcohol (25:24:1) were added, vortexed and heated at 56°C for 15 
min (16). The mixture was centrifuged at 10,000 x g for 10 min and 
the supernatant was transferred into a new tube and processed 
by the silica/guanidinium isothiocyanate method (2). The RNA 
was eluted from the silica pellet with 50 pi of diethyl- 
pyrocarbonate (DEPC) treated sterile water by 15 min incubation 
at 56°C and centrifuged at 10,000 x g for 10 min. The supernatant 
fraction was kept at -20°C until further use. 

Specific oligonucleotide primers that flank the sequence 
corresponding to S protein cleavage site were designed using 
the Gene Runner program version 3.05 (Hastings Software Inc., 
Hastings, NY) (http://www.generunner.com). The primers 
sequences (positions calculated from the start codon of the S 
gene) were: SPK7_F: 5’-TAACTCTTCCGAACCAGCA-3’ (nt 
2085-2103) and SPK7_R: 5'-AATCGCTTCCTAAACAACC-3’ 
(nt 2701-2719), which amplify a predicted fragment of 636 bp. 

The reverse transcription (RT) reaction was performed with 
a mixture of 8 pi of extracted RNA and 2 pi of the random primer 
pdN6 (GE Healthcare, Little Chalfont, UK) and incubated at 
97°C for 4 min. Subsequently, it was placed on ice for 5 min and 
10 pi of RT mix containing lx RTbuffer (50 mM Tris-HCl, pH 8.3, 
3 mM MgCT, 75 mM KC1), 0.1 mM of each dNTP (Invitrogen™ 
Life Technologies, USA), 10 mM DTT, 100 units of M-MLV 
Reverse Transcriptase (Invitrogen™ Life Technologies, USA) 
and ultrapure sterile water to a final volume of 20 pi were added 
and incubated at 37°C for 60 min and followed by inactivation at 
95°C for 5 min. 

For amplification, 8 pi of the RT reaction were added to 42 pi 
of the PCR mix containing 1.5 x PCR buffer (30 mM Tris-HCl pH 


8.4 and 75mM KC1), 2 mM MgCT, 0.2 mM of each dNTP, 1 pi (20 
pmol) of each primer, 2.5 units Platinum Taq DNA polymerase 
(Invitrogen™ Life Technologies, USA) and ultrapure sterile 
water to a final volume of 50 pi. The reaction was performed in 
a thermocycler (PTC-200, MJ Research Co. Water Town, MA, 
USA) with the following time and temperature conditions: one 
step of 4 min/94°C; followed by 40 cycles of 1 min/94°C, lmin/ 
52°C, 1 min/72°C and a final step of 7 min/72°C. 

The products were analyzed by electrophoresis in a 2% 
agarose gel in TBE buffer pH 8.4 (89 mM Tris; 89 mM boric acid; 
2 mM EDTA), stained with ethidium bromide (0.5 pg/ml) and 
visualized under UV light. 

The PCR amplicons were purified using GFX PCR DNA and 
Gel Band Purification (GE Healthcare, Little Chalfont, UK) and 
sequenced in MegaBACE 1000/Automated 96 Capillary DNA 
Sequencer (GE Healthcare), according to the manufacturer’s 
instructions. Sequencing was performed in both directions using 
the forward (F) and reverse (R) primers. Sequences quality 
analysis was performed using Phred/Phrap/Consed Analysis 
Program (http://www.phrap.org) and sequence similarity 
search was performed using BLAST software (http:// 
www.ncbi.nlm.nih.gov/BLAST/). The nucleotide sequences of 
the wild-type Brazilian BCoV strains were aligned and compared 
with the BCoV reference strains using the CLUSTAL/W program. 

The sequences reported in this study have been deposited 
in the GenBank database under accession numbers: 
DQ479421 (BR-UEL1), DQ479422 (BR-UEL2) and DQ479423 
(BR-UEL3). The predicted secondary structure, hydrophilicity 
profile and protease map were determined using Protean of 
DNAStar software (DNAStar, Inc., Madison, WI, USA). 

The amino acid sequence of the cleavage site of the three 
wild-type Brazilian BCoV strains was KRRSSR, indicating an 
amino acid mutation at position 767 (R—> S). Divergence in amino 
acid sequence of the S protein cleavage site among 
coronaviruses strains are indicated in Fig. 1. 

The only mutation described in the proteolytic cleavage site 
of bovine coronaviruses was reported in 1991 with the wild type 
French BCoV strain F-15, which also showed an amino acid 
substitution at position 767 (R—TV) (Fig. 1). However, the authors 
did not discuss the effects related to the amino acid change (24). 

Sequence analysis suggested that the BCoV BR-UEL 
sequences were trustworthy and not the result of RT-PCR or 
sequencing errors, since the same mutation was successfully 
amplified and sequenced in all samples and also in different 
fecal aliquots from each animal. Furthermore, the mutations 
associated with viral adaptation in cell culture were excluded, 
since the samples did not undergo cell culture passage prior to 
PCR amplification and sequencing. 

Although there are no study related to mutational analysis 
in the proteolytic cleavage site of BCoV strains, this 
phenomenon has been previously described for MHV strain. 
Computational analysis demonstrated that the HCoV-OC43, 
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Figure 1 . Predicted amino acid sequences for the spike gene in the cleavage site proximities. Residues in bold represents proteolytic 
cleavage signal sequence of the Mebus strain. The arrow between amino acids 768 and 769 indicates the cleavage site between S1 
and S2. Shaded residue indicates the mutation in the cleavage signal. BCoV BR-UEL strains (GenBank accession numbers 
DQ479421, DQ479422, DQ479423), Mebus (U00735), BCoV ENT (AF391541), LSU94 (AF058943), OK514 (AF058944), Kakegawa 
(DQ479424), LY138 (AF058942), Quebec (AF220295), F15 (D00731), Vaccine (M64668), L9 (M6466). 


MHV-A59 and BCoV presented high degree of identity among 
group 2 coronaviruses (18). Based on these evidences, the MHV 
was used as an appropriate model to hypothesize the effects of 
the new mutation observed in the S protein cleavage signal of 
wild-type BCoV strains during this study. 

Cleavage of S protein is related to efficient cell-cell fusion 
(syncytia) by MHV in infected cell monolayers but is not 
necessary for virus-cell fusion (infectivity). The MHV spike 
proteins, that have mutations that eliminate cleavage into S1 
and S2 subunits, carry out cell-to-cell fusion very inefficiently; 
however, they are able to mediate entry into susceptible cells 
with similar efficiency as wild-type virus (3,8,14). Although the 
cleavage of the S protein is not a prerequisite for fusion, it does 
enhance the induction of syncytia (3). Hingley el al. (12) and 
Hann el al. (10) reported that the cleavage of the S protein of 
MHV-A59 mutant was greatly reduced in comparison with the 
wild-type virus. The sequence at the predicted cleavage site of 
MHV-A59 wild-type is RRAHR while that of MHV-A59 mutant 
is RRADR. The MHV-A59 mutant showed a substitution of a 
weakly basic histidine (H) by negatively charged aspartic acid 
(D) in the cleavage signal (8). These authors suggest that the 
introduction of a negatively charged amino acid into this highly 
basic region destroys the signal or otherwise prevents its 
recognition and thereby inhibits cleavage of the spike. 


Similar to MHV, BCoV S gene also has a proteolytic cleavage 
site formed by a group of basic amino acids upstream from the 
S2 amino terminal region. Except for BCoV F15, this sequence 
is highly conserved in all known BCoV strains (KRRSRR). In 
the BCoV BR-UEL samples was identified an amino acid 
substitution (R—»S) at amino acid position 767 (KRRSSR) (Fig. 
1). This amino acid change occurred due to a single nucleotide 
mutation at the proteolytic cleavage site, CGT to AGT. This 
mutation (C to A), that resulted in the change of the arginine 
codon to serine, had never been described. The cleavage into 
the two subunits occurs between amino acid 768 and 769 and 
is thought to be mediated by cellular trypsin-like proteases 
(19). The protease map indicated that three BR-UEL strains 
lack one predicted cleavage site by trypsin when compared 
with other known BCoV. In MHV studies there are strong 
evidence that furin enzyme is the cell protease responsible for 
cleavage of the S protein in cultured cells (10). In addition, 
although it did not alter the predicted secondary structure or 
hydrophilicity of the protein, the substitution of arginine to 
serine is a nonconservative substitution. Nonconservative 
amino acid changes within a single protein may result in 
alterations of the physical energy of the protein-protein 
interaction and may destabilize the protein native conformation 
(6). Navas-Martin et al. (15) suggest that the proper spatial 
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arrangement of the SI and S2 subunits is crucial for the 
biological functions of the S protein. 

Yoo and Deregt (23) generated BCoV mutants and confirmed 
that a single point mutation at domain II of SI protein was 
responsible for the escape of BCoV from immunological 
response. Similarly, He et al. (11) described that a single amino 
acid substitution in the receptor-binding domain of SARS 
coronavirus S protein disrupted the antigenic structure and 
binding activity. In addition, a single amino acid change has 
already been demonstrated to influence MHV ability to spread 
within the central nervous system (22). 

Recently, Navas-Martin et al. (15) have demonstrated that a 
single cleavage signal substitution in MHV mutants may play a 
major role in the virulence. These authors have associated this 
mutation with higher viral load and highly virulent phenotypes. 
Interestingly, the BR-UEL strains were obtained during an 
outbreak of neonatal diarrhea that culminated with the death of 
calves in a Brazilian dairy cattle herd. Because BCoV infections 
resulting in high mortality is uncommon, the results described 
by Navas-Martin et al. (15) for MHV strain should be also 
investigated for BCoV proteolytic cleavage site mutants. 
Divergence within the S1 gene of others BCoV Brazilian strains 
had already been described by Brandao et al. (4). These authors 
repotted the first description of a gap of 18 nucleotides (deletion 
of 6 amino acids) within the hypervariable region of the SI 
subunit from their isolates and that was also found in human 
coronavirus strain OC43. 

We are uncertain if this single mutation prevented the normal 
proteolytic cleavage of the S protein or increased viral load, as 
occurred in MHV mutants. Therefore, additional experiments 
must be done to determine the effects of the amino acid change 
(R—>S) within the proteolytic cleavage site on the structural 
and functional features of the BCoV S protein. Development of 
a system which introduces infectious cDNA clones with specific 
mutations into the BCoV genome will provide an important tool 
to determine the role of these BCoV mutants in the pathogenesis 
of the neonatal calf diarrhea. 
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RESUMO 

Identificagao de uma mutagao no sitio de clivagem da 

protefna da espfcula em amostras brasileiras de 
coronavirus bovino 

A protefna da espfcula (S), uma glicoprotema de membrana 
do tipo I, e primariamente responsavel pela entrada do virus em 


celulas susceptfveis por meio da interagao inicial com receptores 
celulares especfficos e subseqiiente mediagao da fusao vfrus- 
celula. A protefna S do coronavirus bovino (BCoV) e clivada em 
duas subunidades: a SI, na regiao N-terminal e a S2, na regiao 
C-terminal. O sftio de clivagem proteolftica da protefna S e 
altamente conservado entre as estirpes de BCoV e esta situado 
entre os aminoacidos 763-768 (KRRSRR). Este estudo descreve 
uma mutagao no sftio de clivagem da protefna S de tres estirpes 
do BCoV detectadas em amostras fecais diarreicas de bezerros 
naturalmente infectados no Brasil. O seqiienciamento dos 
produtos de PCR identificou a seqiiencia de aminoacidos 
KRRSSR no sftio de clivagem de nossas amostras, indicando 
uma mutagao na posigao 767 (R—>S). Esta mutagao ocorreu 
devido a uma linica substituigao de nucleotfdeo no sftio de 
clivagem proteolftica, alterando o codon CGT para AGT. Esta e 
a primeira descrigao desta mutagao de nucleotfdeo (C para A), 
que resultou na substituigao do aminoacido arginina por serina 
no sftio de clivagem da protefna S. Neste estudo tambem sao 
sugeridos os provaveis efeitos desta mutagao no sitio de 
clivagem proteolftica utilizando o coronavirus da hepatite dos 
camundongos (MHV) como um modelo comparative. 

Palavras-chave: BCoV, seqiienciamento, protefna da espfcula, 
sitio de clivagem, mutagao 
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